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We provide a unified geometrical description for analyzing the stability of holonomic quantum 
gates in the presence of imprecise driving controls (parametric noise). We consider the situation 
in which these fluctuations do not affect the adiabatic evolution but can reduce the logical gate 
performance. Using the intrinsic geometric properties of the holonomic gates, we show under which 
conditions on noise's correlation time and strength, the fluctuations in the driving field cancel out. In 
this way, we provide theoretical support to previous numerical simulations. We also briefly comment 
on the error due to the mismatch between real and nominal time of the period of the driving fields 
and show that it can be reduced by suitably increasing the adiabatic time. 



PACS numbers: 03.67.Lx, 03.65.Vf 



I. INTRODUCTION 

Recently, there has been a renewed interest for geo- 
metric phases and their application to quantum informa- 
tion [IHtJ including several solid state experiments [HHI . 
These new results could open the way for the realiza- 
tion of "holonomic quantum computation" in which the 
quantum information is manipulated only by means of 
geometric operators. This interesting line of research has 
been opened by the seminal paper of Wilczek and Zee 
[lH and later on by Zanardi and Rasetti [HI, [l4| . 

Holonomic computation is based on an Hamiltonian H 
depending adiabatically and periodically on time via a set 
of time dependent parameters x(t) = (xi(i), . . . ,x n (t)). 
Interesting implementation proposals were initially based 
on atoms and ions [Til [l6| and later in many other quan- 
tum systems |17H23 |. All these models have an Hamilto- 
nian of the form 



H(x) = f(r)H(x) , 



(1) 



where / is a real valued function of r, the norm of a; € M. n , 
and x is the unit vector x/r. This structure in atomic 
physics is usually referred as tripod Hamiltonian. 

The logical operation associated with ([1]) depends only 
on the shadow that the curve x = x(t) projects onto 
the surface of the unit sphere §" _1 (more specifically, 
on the solid angle spanned by it); moreover, the depen- 
dence of the Hamiltonian on f(r) can be transformed 
away by suitable (time-dependent) projective transfor- 
mation in the Hilbcrt space of the system. Thus the 
norm r = r(t) does not need not be periodic — only pe- 
riodicity of x = x(t) is required — and this, in its turn, 
ensures that any mismatch between r(0) and r(T) does 
not affect at all the performance of the gate. 



For sake of concreteness, one may think of the spe- 
cific implementation proposal for quantum dot driven by 
ultrafast lasers (20|, modulated in amplitude and phase 
with the quantum information stored in the excitonic de- 
gree of freedom. However, the same formalism and re- 
sults hold for all the other proposals (l5l422l ]. 

For this purpose, the relevant Hamiltonian is 



H{t) = x{t)-b 



(2) 



with x{t) £ M 3 and suitable matrix- valued vector b. 
More precisely, we have the following structure: 

1. The quantum dot has a level structure of three ex- 
cited degenerate states \i) (i = 1,2,3) at energy 
e, and a ground state |0), set for convenience at 
energy 0. 

2. The system is driven by time-dependent laser fields, 
with frequency in resonance with e, inducing tran- 
sitions between ground and excited states. In the 
interaction representation, the Hamiltonian govern- 
ing the dynamics of the system is ((2J, where b is 
the matrix-valued vector with components 



6 i = |0)<»| + |»)(0|. 



(3) 



The components of the vector x = x(t) repre- 
sent the amplitudes of the three laser driving fields 
which are slowly varying functions of time. More- 
over, the motion of the unit vector x(t) is periodic 
of period T, i.e., 



x(0) =x{T) = x . 



(4) 
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Note that we are not assuming periodicity of r(t). 

We refer to the final geometric transformation as the 
gate operator G, which is defined as the adiabatic limit of 
the evolution generated by ©, for suitable logical space 
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L C C 4 . (We recall that "logical space" stands here for 
"qubit" or, cquivalently, for "two dimensional complex 
space C 2 "), 

G = ad-lim U(T)\ h , (5) 

where U (t) is the evolution operator in C 4 generated by 
([2|), i.e., the solution of (H = 1 from now on) 

.dU 



i-£=H(t)U, U(0)=I, 



(6) 



where / is the identity operator. 

The main advantage in using the gate operator G to 
manipulate quantum information is that it has an intrin- 
sic robustness under different kinds of errors. In fact, 
its performance can be optimized in p resence of environ- 
mental effect and decoherence [23 - [3l| . Similar properties 
and geometric interpretations have been discussed for the 
Abelian geometric phases for which a full analytical treat- 
ment exist s 1321 . These includes environmental induced 
geometric |33ll and non-adiabatic non-Markovian contri- 
butions [3J]. Moreover, numerical investigations show 
that geometric operators are robust against fluctuations 
of the driving fields x(t) [35T4371 ] . so-called parametric 
noise. However, no analytic treatment has been provided 
so far. 

Here we provide a first-principle explanation of such a 
robustness. We analyze the stability of the gate opera- 
tor G under different sources of error which can decrease 
the logical gate performance. We identify two kind of 
errors induced by the fluctuations of the driving fields: 
the error in switching off the driving fields and the one 
accumulated during the evolution. These errors are sup- 
posed to be weak enough to avoid any loss of adiabaticity, 
i.e., they do not produce any transition between non- 
degenerate states, but they can decrease the logical gate 
performance. Using the intrinsic geometric properties of 
the holonomic gates, we show that both errors can be 
reduced in the adiabatic limit, thus recovering the ro- 
bustness of the holonomic transformation. 

The paper is organized as follows. In Section [II] we set 
the stage for the study of the stability of G. In particu- 
lar, we study the adiabatic regime using the language of 
scaling limits and introduce the dimcnsionless adiabatic 
scaling parameter e. In Section Hill we highlight the ge- 
ometrical aspects of the problem, in particular we stress 
the utility of adopting a geometrical extrinsic point of 
view in describing the geometrical features of G. In Sec- 
tion IIVI we discuss the robustness of the holonomic gate 
operator in the presence of parametric noise. In Section 
[V]we conclude. 



II. ADIABATIC REGIME 

In this section we solve (|6]) in three steps. 
First, we observe that H(t) in © is a unitary trans- 
formation of H (0) multiplied by a scale factor, i.e., 



To see how this comes about, consider the time depen- 
dence of the unit vector x = x (t) expressed in terms of 
spherical coordinates 9 = 9(t) and </> = <f>(t) with respect 
to an orthonormal basis i = (ii, i2, 13) in the space R 3 of 
the driving parameters, i.e., 

x(t) = sin 9 cos <j> ii + sin 9 sin <fi 12 + cos 9 '13 , (8) 

and regard R 3 as embedded in C 4 according to the iden- 
tifications ii = |1) , i-2 — |2) , i3 = 1 3) , to which one may 
add, for uniformity of notations, the stipulation io = |0). 
Let 



cos 9{t) cos 4>{t) cos 9 (t) sin <j)(t) — sin9(t) 

— sin^(t) cos </>(£) 

sin 6(t) cos <f>(t) sin 9 (t) sin 4>(t) cos 9(t) 



9{t) 

in the (Ii, 12,13) basis, and 



D(t) = 



1 

<3(t) 



(9) 



(10) 



in the (io, ii, i2, is) basis. Then one may easily check that 
(O is satisfied for 



and scale factor 



R(t) = D{Q)- l D{t) 



a(t) = r(t)/r(0) . 



(11) 



(12) 



The second step consists in writing the equation of 
motion in the moving frame associated with R(t). This 
is realized by the change of variables 



V{t) = R(t)U(t) , 
whence by © and (|7|), 
dV 



dt 



A(t)V + B(t)V, 



(13) 



(14) 



with 



A( t) = i^Rar 



and B(t) = a(t)H(Q) . (15) 



The third step exploits the physical assumption that 
the external fields are slowly varying functions of time. 
Since the energy scale of ([2]) is determined by r to charac- 
terize the evolution we introduce the adiabatic parameter 



1 

Jt 



(16) 



H{t) = a(t)R(t)- 1 H(0)R(t) 



(7) 



where f is an estimate of the size of r(t) (say, the min- 
imum value of r(t) during the drive). The parameter e 
must, of course, be small in order to avoid transitions 
between non-degenerate states. 

For such an adiabatic regime the time dependence of 
the curve in R 3 should be better regarded as given by x = 
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x(et), with the function x(t) having a scale of variation of 
order one. The time for the application of the geometric 
transformation is of order 1/e. Note that the derivative 
of R(t) leads to the rescaling A(t) — > eA(ei), so that the 
dynamical problem in the adiabatic regime becomes that 
of solving the rescaled equations of motions for V(t) 



.dV(t) 
' dt 



eA(et)V(t) + B{et)V(t) 



(17) 



for s <C 1. 

At any given time the Hamiltonian (J2|) can be diago- 
nalized (the explicit time dependence plays no role in the 
diagonalization and we shall omit it in the notations). 
Consider a vector u in C written as u = uq\q + u, where 
u = mil + u 2 i2 + u 3 i 3 and x = xi'n + X2I2 + X3I3. 
From the definitions of ij and bj in we obtain 
bjik = (Jj/Jo + Sokij- Immediately from the eigen- 
values equation reads 



H [u i + u] = x ■ u i + u x , 



(18) 



from which eigenvalues and the eigenvectors of H follow: 
take u = x, then from 03 



H[u i +x]= r[i + u Q x] . 

Thus, Uo = 1 gives the eigenvalue A + = r and uq = 
— 1 gives the eigenvalue A_ = — r. The corresponding 
normalized eigenvectors are, respectively, 



e± = (1/V2)(x±i ). 



(19) 



Finally, for uq = and u orthogonal to x (as vectors in 
M 3 ), we read from (fTBf that Ao = is a doubly degenerate 
eigenvalue. Thus, we have recovered well known proper- 
ties of the tripod Hamiltonian @ : there are two degener- 
ate states, at zero energy, called "dark states" and other 
two, called "bright states" , with one excited state (e + in 
(Unj)) and one ground state (e_ in p^|) ) with energy r 
and — r, respectively. In the following, we shall denote 
by Pp, fi = +1,-1,0, the spectral projectors of H(0) 
corresponding respectively to the eigenvalues A_|_ = r(0), 
A_ = — r(0) and Ao = 0. Po projects onto the plane L 
orthogonal to x(Q) which soon will be identified with the 
logical space in (|5|). 

Consider now the interaction representation of the time 
evolution operator V(t), whose dynamics is governed by 
(|17p , with respect to free dynamics generated by B 



Vi(t) = W(t)V(t), 



(20) 



with 



W(t) = e'fo B ^ dt ' = e ^*(0) = J2 e 4 ^^ Pp , 

P 

(21) 

where (sec 02])) 



h(t) = / a{r)dT . 
'0 



(22) 



Then from Q7|) it follows that V/(i) satisfies 



i^-jp- = eWWAWWityiVitf) , (23) 



or, equivalently 



Vi(t) =1+1 W(s/s)A(s)W(s/e)- L V I (s/e)ds. (24) 
Jo 



Equation ([24]) shows that the effect of A on the evo- 
lution manifests itself only on the adiabatic time scale 
t e = t/e. Thus, setting Vf(t) = Vi(t/e) one gets 

pet 

Vf(t)=I + J2 / e i ° {Xf> ~ X ' , ' ) Pi}A{s)Pp,V?{s)ds. 



, Jo 



(25) 

Since h(t) is positive and never equal to zero, standard 
stationary phase approximation gives /?'=/? and 

Vf{t) =I + J2 [ PpA(s)PpV I e (s)ds + 0(e) . (26) 

a JO 



By multiplying both sides of the above equation by P a we 
see that the evolution separates into autonomous spectral 
components. In particular, since Po projects onto the 
plane L, (|26[) defines a non-trivial dynamics on it. In 
other words, the operator 



G(t) = P VKt)P 



(27) 



evolves autonomously in the adiabatic limit (modulo cor- 
rections of order e) according to the equation 



where, recalling p"5]) . 



A(t) = iP ^-R(t)- 1 P . 



(28) 



(29) 



III. GEOMETRY 

Equation (]2^)) can be solved analytically using a geo- 
metric approach. By defining the operator-valued vec- 
tors a = (01,02,03) and a* = (0^,02,03), having com- 
ponents Oj = \0)(i\ and a* = \i}(0\, i = 1,2,3, re- 
spectively. The Hamiltonian in ([5]) can be written as 
H(t) = a ■ x(t) + a* ■ x(t). Another operator-valued vec- 
tor that it is useful to introduce is J = a* x a, whose 
components 



Ji 

J-> 



a 2 a 3 ~ a 3 a 2 

a^ai — 0*03 
0^02 — O^Oi 



(30) 
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are indeed the generators of an SO(3) algebra with com- 
mutation relations [Ji,Jj] = —eijkJk and in terms of 
which D(t) in (fT0|) can be expressed as 



D(t) 



- P -e(t)J2 p -<P(t)J3 



(31) 



The frame e(t) = (e , eg(i), e^i), e r (i)), where 



eo = io 

e r (t) = x(t) 

eg(t) = cos 8(t) cos (f>(t) ii + cos 6(t) sin cf>(t) '12 — sm9(t) 13 

e 4>{t) = — sin0(t) ii + cos0(t) i2 (32) 

is indeed a moving frame adapted to the surface of the 
unit sphere § 2 on which x(t) moves in the course of time, 
i.e., eg(t) and e<p(t) are tangent to S 2 , and e r (t) is per- 
pendicular to it. Then pip is the operator transforming 
the frame i = (io, ij., '12, 13) into the frame e(t), 

D(t)-\=e e (t), D{t)- 1 i 2 =e ip {t), D(t)-% = e r (t) . 

(33) 

Note that D(0) 1 is not the identity and, in particular, 
that L»(0) _1 i 3 = x(0), whence from ((TTJ) 

x{t) = e r {t) = D(t)- l D{Q)x{Q) = R{t)- l x{Q) . 

Then ([7]) is nothing but an expression of the usual du- 
ality between action of the operators on vectors and on 
operators, i.e., 

x(t) • b = [mt)- 1 ^)] ■ b = sfe(O) • [ii(i)bi?(i) _1 ] , 

andA(0)-b = L>(0)- 1 fe 3 £'(0). 

The unit vectors eg(t) and e^(i) (the "dark states") 
are a natural basis in the moving degenerate space — the 
plane orthogonal to e r (t) = x(t), which is, according 
to Section [ill the degenerate eigenspace of the eigen- 
value of H(t). Let Po(t) denote the projector onto 
such a plane, then Po(0) is the projector P onto the L 
plane. Note that Po is just a D(0)-rotation of the pro- 
jector P(iii2) = a\a\ + a^a^ onto the ii-i2 plane; thus, 
Po = P(0)- 1 P(i 1 i 2 )£»(0). 

Putting all the pieces together, from (fTTj) . ([25)1 and 
(|29]l we obtain 

A{t) = i ^(0)- 1 P(i 1 i 2 )^^P(<)- 1 P(i 1 i 2 )P(0) 

at 



with (see pi 
dD(t) 



dt 



D(t)" 1 = -8J 2 - <j) (sin 6J1 + cos 9J 3 ) 



Using (f3Tj) and the above definition of P(iii2), we can 
calculate the projection P(i!i 2 ) J^P(iii2) = S^Jk and 
we are left with 

A{t) = i{^cos9)D{0)~ 1 J 3 D{0) . 

Observing p"Tj) at time t = 0, and using the rotational 
properties of Jk we have 



and, finally, we arrive to a geometric expression 
A{t) = i(<j)cos9)J ■ x(0) . 



(35) 



We have at our disposal all the ingredients to deter- 
mine the geometric operator that is used to manipulate 
the quantum state. Now, since A(t) at different times 
commutes, equation (J2FJ can be solved by direct expo- 
nentiation 



G(t) = cxp 



1 1 cos 9(t')cj)(t')dt' J ■ x(0) 




(36) 



One may recognize that G(t) is a rotation in the L plane 
of an angle given by the integral multiplying J ■ x(0) 
in p6[) . At this point starting from ([27} and going back 
from p6p to the moving frame by means of equation ([2"0|) 
is immediate: Wit) is just the identity on L since the 
corresponding eigenvalues Ao is zero. To go back to the 
laboratory frame, we follow (|13p and obtain 



ad-lim U(t)\ h = P(0 _1 G(t) . 



(37) 



From the assumption ([4]) of periodicity of x(t), it fol- 
lows that R(T) = I. Therefore 



G = ad-lim U(T)\ L = G(T) . 



(38) 



If we choose a logical basis directly in L we can write 
G(t) in ([36]) as a 2 x 2 matrix. Thus (modulo change of 
basis in L), 



where 



G = G(T) 



n = 



cos il sin Q, 
— sin f2 cos fi 



cos9(t)4>{t)dt. 



(39) 



(40) 



£>(0) _1 J 3 £>(0) = J ■ sc(0) 



(34) 



This is the expected result: the G(T) operator depends 
only on the solid angle spanned by x(t) on the Hamilto- 
nian parameter space. Note that this result does not 
rely on the assumption of full periodicity (x(0) = x(T)) 
made by Wilczck-Zee and Zanardi-Rasetti, although it 
holds only for a restricted class of Hamiltonians. 



IV. STABILITY 

We shall now consider the effects of external perturba- 
tions on the system, in order to understand the robust- 
ness of quantum evolution in presence of noise. Several 
are the sources of undesired energy exchange due to cou- 
pling with impurities, and/or external environments, or 
due to imprecise control of the system parameters dur- 
ing the evolution. In the following we will focus on this 
parametric error. First of all, we observe that the good 
performance of the gate G relics on the validity of the 
adiabatic limit, and since the dimcnsionless constant e is 
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finite, this fact alone introduces an error of order e, i.e., 
0(e), which is the order of magnitude of the off-diagonal 
terms (/3 ^ j3') in (|25jl which are neglected in the limit 
e — > 0. Then we shall say that the gate G is stable if 
the perturbations produce corrections of higher order in 
e. Letting aside corrections not following a power law, 
we can say that stability is ensured if the corrections on 
G due to the perturbations are 0(e r ), with r > 1. To 
simplify the analysis, we will work with dimensionlcss 
quantities: energy scale in units of f, so that (fTB|) be- 
comes e = 1/T. 

We consider the error induced by the inaccuracy of 
the control field. We start with the case in which the 
actual curve in the parameter space is not x(t) but in- 
stead x'(t) — x(t) + 8x(t), fulfilling still the periodicity 
requirement (QJ on the unit vector x'(0) = x'(T). The 
error 5x{t) is assumed small with respect to x(t) in the 
sense of some suitable functional norm. It should be re- 
garded as rapidly fluctuating random process whose scale 
of variation is very small on the adiabatic scale. 

From the above periodicity it follows that the holo- 
nomic operator U'(T), related to x'(t) is still of the form 

ad-lim U'(T)\ h = G\T) = G[x'(t)] . (41) 

The geometric operator G[x'(t)} is now a functional of 
x'(t) and then the error is accumulated during the whole 
evolution. 

In the following we will evaluate at lowest order in 
8x(t) the error on the holonomic operator U'(T) = 
G[x(t) + 8x(t)] induced by the fluctuations of the driving 
fields along the path. We have 

U'(T) = G[x(t) + 8x{t)] = G[x(t)] +8G + o(a) (42) 

where a is a measure of the size of the (mean) variation 
of 8x(t). To do that, let us start to rewrite the solid angle 
in ((301) for path x'(t) as 




where 7 is the shadow on the unit sphere § 2 of the curve 
7 in R 3 given by the parametric equations x' = x'(t), 
and satisfying the conditions ((4]) of partial periodicity, E 
is the surface on § 2 bounded by 7, i.e., 7 = SE. 

Thus, ft is the area of E, that is, the solid angle 
spanned by curve 7 (i.e., the solid angle from which 7 is 
seen from the origin in K 3 ). Accordingly, different uni- 
tary geometric transformations and then quantum logical 
gates can be constructed traversing different loop in the 
parameters space. 

Since the North Pole = <f> = is a singularity of 
spherical coordinates, one is naturally lead to consider 
the 1-form on 8 2 , locally defined by lj = cos 8d(j) and 
extended it to all § 2 in a coordinate independent way. 
Accordingly (|43)) should be replaced by the Stokes the- 
orem on § 2 expressed in an "intrinsic" geometrical way 
(i.e., coordinate independent), 

n = / uj = / duj . (44) 
Jot. Jj: 



This is the a ppr oach usually adopted in holonomic com- 
putation (l3lll4| . To take care of the singularity problem 
one can rewrite (|4"3"]l in terms of the auxiliary vector field 
A = e^(l — cos 9) /(r sin 9) 

Cl[x'(t)] = J A-dr = JJ B-dS, (45) 

where 

B = V x A = ie r . (46) 
We can now rewrite (|45[) as an integral over time 

Q[x'(t)]= [ A(x'(t))-x'(t)dt. (47) 
Jo 

Noticing that the right-hand side is analogous to the La- 
grangian of a particle in a magnetic field given by (|46[) . 
we have (see and (|3"9")l ) 

8G = i[G J ■ x(0)}8Q , (48) 

with 

8Q= [ B(x (t)) X x(t) ■ Sx(t) dt (49) 

Note, as expected, that radial fluctuations give no con- 
tribution to the variation since B x x is tangent to the 
sphere. 

We now define the statistical properties of 8x(t) which 
describes the parametric noise perturbing the external 
field. As already discussed it should have a scale of vari- 
ation very small on the adiabatic scale, yet, if we wish 
to ensure the validity of the adiabatic approximation, we 
should demand, at the same time, that its scale of vari- 
ation be sufficiently long on the microscopic scale. The 
simplest possibility to ensure this is to regard the com- 
ponents of 8x(t) as independent mean-zero stationary 
Gaussian processes with "white-noise" correlation func- 
tion, i.e., 

<8x i (t)8x j (t') >= 6,,r,ajMI - t') , (50) 

with i,j = 1,2,3 and where of is a measure of the 
strength of the noise components (time- independent, as 
the processes are stationary), r^s are the correlation times 
of the noise components. Since the Gaussian processes 
are stationary, we can take of =<<5.Ti(0) 2 >. 

Again, it is convenient to consider n and <7j as func- 
tions of e, e.g., 

n = 0{? ) , tn = 0(e») (51) 

with suitable exponents p > and q > 0. Recalling 
that the intrinsic error of the adiabatic approximation is 
0(e), we should then inquire whether there is a range of 
p, q- values for which the mean error 

a = V<sn 2 > - <sn> 2 (52) 
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is below the O(e) upper bound, i.e., 0(e r ) with r > 1. Of 
course, A would then provide an estimate of the (mean) 
first order correction <5G in (|42j) . 

Since the fluctuations dxi have zero mean, <6Q>= 0. 
Moreover, from (|49l) and (|50|) it follows 



(53) 



3 .T 

< Sn 2 >= V r it 7? / [B x x] 2 dt . 



Recalling (|46|) . we find 



A 



3 -T 
V^of / i 2 

.•— 1 Jo 



(54) 



Since the velocity ii scales as 1/T, the integral in the 
right-hand side of is 0(T x 1/T 2 ) = 0(1/T), i.e, 
0(e), and the mean error is then 



A = 0(e r ), r=p/2 + q + l/2. 



(55) 



Prom (fSTj) it follows that with this constraint there are 
many solutions p/2 + q+ 1/2 > 1, in the desired range of 
p, g-values, which ensure stability of the gate. Moreover, 
from ([55)1 , we also read the answer to the questions about 
the smallness of and largeness of oi , namely, that the 
fluctuations can indeed be quite large with respect to e, 
provided that the t^s are small on the macroscopic scale 
(but large on the microscopic scale). This is the can- 
cellation effect already discussed in Rcfs. [35|, [36|. For 
example, let e = 10~ 4 . Then a noise with correlation 
time of order, say, t = 10 -2 (p = 1/2) can have fluctua- 
tions of order, say, a = 10~ 2 (q = 1/2) while producing, 
at the same time, an error on the gate of order 10~ 5 
(p/2 + q + 1/2 = 5/4), well within the bound e = 1(T 4 
of the adiabatic regime. 

We conclude this section with two observations. The 
first concerns the geometric interpretation of solid angle 
perturbation, with a simple geometrical formula for the 
right-hand side of ([5"3"|) . We focus on the case of noise 
components 5xi(t) with the same statistical properties, 
i.e., Ui = a and Tj = r (a condition which is indeed quite 
reasonable from a physical point of view) . Since is in- 
variant under re-parametrization of time, it is convenient 
to use as invariant parameter the arc length s with origin 
in Xq and rewrite (|4"Tf as 



fl[x'{s)} = / A{x'(s))-x'{s)ds, 
Jo 

where L is the length of the curve and x = dx/ds. Then 
([53)) becomes 

3 ,L 

< 5n 2 >= la 2 V / [B x x] 2 ds , 

where I is the correlation length of the noise. (Note that 
the integral is now 0(1), since L is the length of the 



curve on § 2 , and therefore i = 0(e p+1 )). But, from (|46|) 
B = n, the normal to the curve lying on S 2 , i.e., 



<Sn A >= £er 



|n x x\ds , 



and one recognizes 



SA = a |n x x\ds 
Jo 

as the area of a thick boundary along the curve of width 
a. Thus, A 2 = lad A ■ 

Let us now turn to the second observation. It con- 
cerns how to deal with the difference between nominal 
and real value of the period of the laser. We recall that 
in the tripod model two control fields are turned off at 
the initial and final configurations [l5Tj22j . However, in a 
practical implementation, there is uncertainty in the real 
time T in which the fields are turned off since in general 
this does not correspond to the nominal time To, with 
T = T + AT, being AT a statistical fluctuation time. 
This can be due, for example, to latent times or delays 
in the control fields of the experimental set-up. In the 
parameter space this error corresponds to an evolution 
that, if referred to the nominal time, has i? _1 (To) 7^ /, 
since the periodicity requirement is fulfilled at time T, 
e.g., i? _1 (T) = I. This means a path non-periodic at the 
nominal time (391. Thus, from (|37| 



(56) 



f/(T ) = i?- 1 (T )G(T ), 



and therefore at the lowest order in AT 



R-^Tq) = iT x (T - AT) w / 



dR- 



dT 



AT. 



(57) 



For the partial derivative in the right-hand side of the 
previous equation, one has 



dR- 



dT 



dR- 



<)x 



x(T). 



(58) 



The <9i? -1 /dx depends only on geometrical properties of 
the path and thus is of order one, while the velocity is 
of order 1 /T (since we are analyzing the system in the 
adiabatic regime). Thus the error due to the mismatch 
between real and nominal time of the period of the laser 
is of order 1/T w 1/T) and can therefore be reduced by 
suitably increasing the nominal time Tq. 



V. CONCLUSIONS 

We have provided a unified geometrical description for 
analyzing the stability of holonomic quantum gates in 
the presence of parametric noise affecting their time evo- 
lution. We have identified two main critical parameters: 
the correlation time of the noise and its strength with 



7 



respect to the driving field. In this way, we have recov- 
ered what was already obtained numerically for Abelian 
(35l ] and non- Abelian gates [IB] , namely, that even strong 
fluctuations in the driving field can lead to accurate logi- 
cal gates if the correlation time is short enough to let the 
fluctuations cancel out. This is an effect of the geometric 
dependence of the holonomic operator. In addition, we 
have shown that the error due to the mismatch between 
real and nominal time of the period of the laser can be 
reduced by suitably increasing the adiabatic time. 
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